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PREFACE 

This  report  is  one  in  a  series  of  AEDC  Technical  Reports  which  document  the  solu¬ 
tions  of  the  Orr-Sommerfeld  and  Rayleigh  equations  and  how  these  fluctuations  are 
generated.  The  reports  in  this  series  are: 

(1)  Rotational  and  Irrotational  Freestream  Disturbances 
Interacting  Inviscidly  with  a  Semi-Infinite  Plate 
(AEDC-TR-83-3) 

(Z)  Exponentially-Varying,  Standing  Waves 
in  Parallel-Flow  Boundary  Layers 
(AEDC-TR-83-4) 

(3)  Waves  Which  Travel  Upstream  in  Boundary  Layers 

(AEDC-TR-83-7) 

(4)  Spatially-Decaying  Arrays  of  Rectangular  Vortices 
Interacting  with  Falkner-Skan  Boundary  Layers 

(AEDC-TR-83-8) 

(5)  The  Boundary-Value  Problem 

for  Two-Dimensional  Fluctuations  in  Boundary  Layers 
(AEDC-TR-83-9) 

(6)  Nonperiodic  Fluctuations  Induced  by  Stationary  Surface  Waviness 
on  a  Semi-Infinite  Plate  (this  report) 

(AEDC-TR-83-10) 

Reports  (5)  and  (6)  are  concerned  with  unsteady  freestream  disturbances  and  with 
stationary  surface  waviness,  respectively.  However,  the  two  studies  are  closely  related 
mathematically.  If  the  solutions  are  properly  expressed,  then  the  solution  for  (6)  is 
imbedded  within  the  solution  for  (5)  when  the  frequency  vanishes  and  when  the  surface 
waviness  is  in-phase  on  both  sides  of  the  plate.  Hence,  an  analogy  exists  between  these 
physically  quite-different  cases. 

All  reproducibles  used  in  the  reproduction  of  this  report  were  supplied  by  the  authors. 
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1.0  INTRODUCTION  AND  LITERATURE  SURVEY 
1.1  Introduction  and  Motivation  for  this  study 

Aerodynamic  surfaces  are  wavy  in  various  degrees  as  a  result  of  machining 
and  fabrication  of  the  surface.  The  deformation  of  the  structure  due  to 
aerodynamic  loading  and  body  forces  can  also  cause  wrinkling  and  local 
buckling  of  the  skin.  Subscaled  models  built  for  aerodynamic  testing  are 
particularly  subject  to  adverse  affects  of  surface  waviness,  even  when  the 
aerodynamic  surfaces  are  manufactured  with  very  small  tolerances,  because  the 
waviness  can  be  large  nondimens ional ly.  Unsteady  surface  waves  may  arise 
through  the  use  of  active  or  passive  compliant  surfaces,  panel  flutter,  a 
layer  of  water  flowing  along  the  surface  having  a  wavy  air/water  interface, 
etc.  While  the  theory  we  develop  in  this  study  can  be  easily  adapted  to  the 
case  of  traveling  waves,  the  problem  we  consider  has  stationary  waviness. 

The  surface  waviness  can  affect  laminar  and  turbnlent  skin  friction, 
transition  to  turbulent  flow,  and  the  structure  of  the  boundary  layers.  If 
the  pressure  along  the  surface  is  not  the  same  at  equal  heights  of  surface 
displacement,  h(x) ,  then  the  waviness  will  cause  pressure  drag.  The  waviness 
(on  an  otherwise  flat  plate)  can  lead  to  local  separations  after  only  a  few 
wavelengths.  At  the  very  high  Reynolds  number  flows  of  aerospace 
applications,  Gortler  instabilities  are  possible  in  the  regions  of  streamline 
concavity. 

While  it  is  nseful  to  analytically  assume  that  the  wavy  wall  extends  to 
infinity  upstream  and  dovnstream,  the  wavy  region  may  end  at  leading  and 
trailing  edges,  or  the  wavy  wall  may  connect  with  (or  flare  into)  smooth 
walls.  While  it  is  useful  to  numerically  study  flows  which  are  spatially 
periodic,  the  inflow  and  outflow  boundary  conditions  will  not  be  the  same  if 
edges  or  flaring  regions  are  "nearby".  The  nonperiodicity  of  the 
inflow/ontf low  boundary  conditions  will  influence  the  flowfield,  even  though 
the  surface  displacement  may  be  periodic.  Even  the  periodic  inf low/outf low 
conditions  can  be  altered  by  edges  of  the  wavy  regions. 

The  case  of  inviscid,  irrotational  flow  over  a  small- amplitude  sinusoidal 
surface  was  analyzed  by  Kelvin  and  Helmholtz  (Ref.  1}.  In  a  coordinate  system 
where  the  surface  waves  are  steady,  and  introducing  a  phase  shift  0  which  is 
useful  later,  the  classical  solutions  for  the  longitudinal  and  normal  velocity 
fluctuations  are 

u(x,y)=u<)SiIi(ax-0)exp(-ay) 

v(x,y)=v#cos{ax-9)exp(-ay)  (l.la.b) 

However,  generally  i£  leading  or  trailing  edges  are  present .  the 
Kelyin-HelTrbp) tm  solutions  are  not  mathematically  complete :  the  flow  cannot  be 
represented  as  a  superposition  of  the  Kelvin-Helmholtz  solutions.  Other 
waveforms  may  be  required  in  the  analysis.  Even  for  the  steady  case  of  an 
inviscid  uniform  mean  flow,  a  mathematically  complete  set  includes  a  pair  of 
steady,  exponentially-varying  standing  waves  of  form 
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u(x,y)=cos(0y)exp(±0x) 

v(x,y)=sin(py)exp(±0x)  (1.2a,b) 

where  0  is  a  real  number.  These  waves  grow  or  decay  exponentially  in  the 
streamwise  direction  and  satisfy  impermeability  along  a  flat  plate.  While 
there  also  is  the  possibility  of  another  family  of  waves  representing  vortical 
fluctuations  in  the  freestream,  and  some  other  variations  to  (1.2a,b),  the 
above  forms  are  sufficient  for  this  discussion.  We  are  particularly 
interested  in  the  decaying  form  of  these  standing  waves  and  how  they  are 
linked  to  the  Kelvin-Helmholtr  solutions  through  the  presence  of  the  leading 
edge.  The  issue  of  mathematical  completeness  is  further  discussed  in  Section 
1.6. 


Solutions  (1.2a,b)  are  simplified  versions  of  waves  of  form 
u(x,y)=f(y)exp(±px— iut) 

v(x,y)=0(y>exp(±px-iu>t)  <1.3a,b) 

found  in  viscous,  parallel-flow  boundary  layers.  These  waves  have  been 
calculated  by  Rogler  and  Tsugd  (Ref.  2),  and  Rogler  (Refs.  3a, b)  as  solutions 
of  the  Orr— Somraerf eld  equation  for  Falkner-Skan  boundary  layers. 

A  motivation  for  this  study  is  to  better  understand  the  role  and 
properties  of  these  additional  solutions.  Even  for  cases  where  superpositions 
of  the  Ke 1 vin— Helmhol tz  solutions  are  possible  (e.g.  the  case  of  a 
semi-infinite  flat  plate  joined  to  a  semi- inf inite  sinusoidal  wall)  it  may  he 
useful  to  represent  the  flowfield  in  terms  of  mathematically  complete  sets  in 
the  two  regions,  and  join  them  smoothly  at  the  junction.  An  ob iective  of  this 
study  is  to  obtain  the  flowfield  over  a  semi-infinite  sinusoidal  wall.  While 
we  consider  the  wavy  surface  to  be  stationary,  with  appropriate 
uondimensional izations  and  phase  shifts,  the  solution  for  the  flowfield  is 
also  applicable  for  flow  over  a  semi-infinite  plate  with  traveling  surface 
waviness.  The  theory  with  traveling  surface  waviness  is  documented  in  Ref .22. 


1 .2  An  introduction  to  the  influence s  of  leading  and  trailing  edges  of  wavv 
surfaces  and  flaring  sections  between  wavv  and  smooth  surfaces 

Several  investigators  have  been  aware  of  (a)  the  influence  of  the  phase 
angle  0  when  a  leading  edge  is  present,  (b)  the  length  of  the  flat  section 
which  joins  to  the  wavy  region,  and  (c)  geometrical  details  of  the  flaring 
section.  There  are  several  mechanisms  by  which  the  leading/trail ing  edges  or 
flaring  sections  can  influence  the  flowfield.  For  simplicity,  we  will  confine 
this  discussion  to  first-order  boundary  layer  effects,  where  these  factors  can 
influence  the  freestream  velocity  at  the  boundary  layer  edge,  ^(x)  . 

(1)  In  the  doubly- inf ini te  Kelvin— Helmholtz  solution,  the  phase  9  of  the 
surface  waviness  can  be  eliminated  by  a  simple  coordinate  shift.  However,  the 
calculations  by  Dryden  (Refs.  4,6,  based  on  Ref.  5)  for  the  boundary  layer 
development  along  a  plate  show  that  the  separation  point  depends  on  the  phase 
9  of  the  sinusoidally  imposed  freestream  velocity 
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UM(x)  =  U,#  +  esin(ax-6)  (1.4) 

(We  have  freely  changed  the  definition  and  sign  of  the  phase  angle  used  by 
other  investigators  to  be  consistent  with  the  present  wort.)  With  a  leading 
edge  presenti  the  phase  angle  of  the  waviness  becomes  a  relevant  parameter 
because  the  phase  influences  the  initial  condition,  the  imposed  pressure 
gradient  associated  with  the  velocity  (1.4),  and  the  edge  boundary  condition 
in  the  boundary  layer  equations.  Dryden's  calculations  (Ref.  4) 
systematically  considered  these  effects. 

(2)  The  flow  around  the  leading  edge  can  substantially  influence  the  initial 
conditions,  the  pressure  gradient,  and  the  edge  boundary  conditions.  This 
flow  depends  on  the  phase  of  the  waviness.  Based  on  the  theory  presented 
herein,  the  freestream  velocity  takes  the  form 

Dro(x)  =  Uro  +  esin(ax-Q)  +  eu^s^(x,6)  (1.5) 

\ 

' —  Secondary  flow  caused  by  the 
leading  and  trailing  edges, 
flaring  sections,  etc. 

—  Kelvin-Helmholtz  solution 
evaluated  along  y  =  0 

—  uniform  mean  flow  which  would  exist 
if  the  plate  were  flat 

This  expression  is  valid  only  along  y=0  where  a  sinusoidal  wall  is  present. 
The  secondary  flow  is  also  proportional  to  the  amplitude  of  the  surface 
waviness.  The  secondary  flow  is  not  period ic  in  x.  We  shall  show  that  it  is 
composed  of  a  superposition  of  decaying  waves  of  form  (1.2a,b).  While 
would  vanish  far-downstream  of  the  leading  edge,  the  influence  in  the  boundary 
layer  could  persist  downstream  because  of  the  development  of  the  boundary 
layer  in  the  first  half-wavelength  or  so  of  the  waviness.  A  similar  behavior 
is  possible  with  a  flat  section  connected  to  a  wavy  section.  In  elliptic 
problems,  the  (possibly  different)  phases  of  the  waviness  on  the  top  side  of 
the  plate  end  the  bottom  side  of  the  plate  would  he  relevant.  In  the  present 
work,  we  assume  that  the  phases  are  the  same. 

(3)  If  both  leading  and  trailing  edges  are  present .  then  the  waviness  can 
influence  the  circulation  about  the  plate.  The  secondary  flow  will  now  depend 
on  the  plate  length.  L.  This  secondary  flow  in  eqn.{1.5)  could  remain 
significant  for  many  wavelengths  from  the  leading/trailing  edges.  A  similar 
effect  could  arise  with  a  finite-length  of  surface  waviness  along  a 
douhly-inf inite  plate.  The  feature  shared  by  both  of  these  cases  is  that  the 
end  regions  of  the  wavy  surface  can  establish  high  and  low  pressures,  and  the 
difference  between  those  pressures  divided  by  the  length  of  the  wavy  section 
is  a  pressure  gradient  which  could  have  a  relatively  long  range  influence.  In 
elliptic  cases,  another  feature  of  relevance  is  whether  the  model  support  or 
attachment  of  the  wavy  surface  to  the  wind  tunnel  wall  will  influence  the 
circulation  about  the  plate. 

(4)  The  geometrical  details  of  any  (non-sinusoidal )  flaring  sections  which 
connect  the  wavy  surface  to  adjacent  non-wavy  surfaces  could  he  treated  as 
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still  another  superposed  flowfield.  This  flow  would  he  excited  by  the  wall 
boundary  condition  which  is  the  difference  between  the  sinusoidal  wall  and  the 
flaring  region. 

Analytical  solutions  have  been  obtained  (Refs.  7-9)  for  the  interaction 
of  various  freestream  disturbances  with  semi-infinite  plates.  Within  the 
framework  of  inviscid  linearized  airfoil  theory,  the  closely  related  solution 
for  flow  over  stationary  surface  waviness  on  a  semi-infinite  plate  can  be 
developed.  For  this  reason,  a  semi-inf inite  wavy  wall  is  used  in  this 
analysis.  The  closed-form  solution  which  results  can  be  rearranged  to 
illustrate  the  origin  of  exponentially— varying  standing  waves.  In  this 
report,  we  will  not  directly  consider  either  the  combined  effects  of  leading 
and  trailing  edges,  or  the  effects  of  nonsinnsoidal  flaring  regions. 

In  Sections  1.3-1. 9,  we  outline  eight  approaches  to  the  problem  of 
studying  flow  over  wavy  walls.  While  the  survey  is  not  exhaustive,  additional 
references  in  any  category  can  be  found  through  the  papers  and  reports  which 
we  cite.  We  did  not  thoroughly  survey  the  related  topics  of  (1)  flows  over 
traveling  wavy  walls,  (2)  geophysical  studies  of  the  excitation  of  water  waves 
by  the  wind,  (3)  airflow  over  a  flat  plate  with  a  film  of  fluid  with  a  wavy 
surface,  (4)  Gortler  instabilities  in  the  concave  streamlines  over  a  wavy 
surface,  (5)  influence  of  waviness  on  the  Tollmien-Schlichting  instability, 

(6)  the  triple-deck  analyses  of  flows  over  a  small  hump  in  a  boundary  layer, 

(7)  a  layer  of  liquid  flowing  over  a  wavy  surface,  such  as  the  flow  of  a  river 
over  an  undulating  riverbed,  (8)  heat  transfer  in  corrugated  pipes,  (9)  low 
Reynolds  number  lubrication  problems  with  wavy  surfaces,  (10)  flow  in  hlood 
vessels  and  pipes  with  wavy  walls,  (11)  flows  along  flat  eels  and  other 
animals  with  undulating  surfaces,  and  (12)  boundary  layer  development  with 
surface  roughness. 


1 .3  If rotations  1  flows  past  smal 1— amp  1 itnde  and  large— amp 1 itude  wavy  walls 

Irrotational  incompressible  flow  over  a  small— ampl itude  sinusoidal  wall 
was  analyzed  by  Kelvin  and  Helmholtz  (Ref.  1).  The  streamlines  of  this  flow 
are  plotted  in  Figure  1.  This  solution  was  extended  by  Ackeret  (Refs.  10-13) 
to  the  cases  of  compressible  subsonic  and  supersonic  flow.  For  the  subsonic 
cases,  the  pressure  is  symmetric  on  the  windward  and  leeward  sides  of  the 
wave,  and  there  is  no  pressure  drag.  For  the  supersonic  case,  wave  drag  is 
generated  because  of  the  asymmetric  variation  in  pressure. 

Transonic  flow  over  a  wavy  wall  was  analyzed  by  C,  Kaplan  (Ref,  14)  and 
Hosekawa  (Ref.  15)  who  included  the  nonlinear  transonic  term  in  their 
analyses.  Both  assumed  the  flow  to  be  isentropic  and  periodic,  and  thus 
tacitly  neglect  the  increases  in  entropy  across  the  shocks.  Moore  and  Gibson 
(Ref.  16)  and  Vincenti  (Refs.  17,18)  analyzed  inviscid  compressible  flow  over 
a  sinnsoidal  wall  of  a  gas  in  vibrational  or  chemical  non-equilibrium.  They 
focused  on  the  effects  of  finite-rate  chemical/vibretional  adjustments,  and 
neglected  the  effects  of  any  (nonlinear)  transonic  term  in  their  consideration 
of  subsonic,  transonic,  and  supersonic  flows.  Kyhming  (Ref.  19)  discussed  the 
combined  problem  of  transonic  effects  and  finite-rates  of  adjustment.  For  the 
(equilibrium)  transonic  cases,  the  mixed  supersonic/subsonic  flows  cause 
pressure  drag.  For  the  non-equilibrium  cases,  the  finite  rates  of  adjustment 
will  cause  pressure  drag  even  In  the  subsonic  cases. 
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Figure  1.  Streamlines  for  the  periodic  Kelv in-Helmholtz  solution,  y,^th) ^  for 
incompressible,  irrotational  flow  over  a  wavy  wall.  The  phase  of 
the  waviness  in  this  figure  is  the  same  as  in  Figures  2b,  5b, and  6f. 


Figure  2.  Geometry  and  streamlines  for  nonperiodic  flow  in  a  corner  induced  by 
a  traveling  wavy  wall  (Ref.  23). 
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Horlock  (Ref.  20)  analyzed  incompressible  flow  in  an  unsteady  wind  tunnel 
with  small-amplitude  wavy  walls  with  various  phase  relationships  between  the 
two  wavy  walls.  Rogler  (Ref.  21)  analyzed  the  small-amplitude  cases  of  wall 

amplitudes  varying  as 

h(x,t)  =  h0es;p[ia(i-ct)]  with  a  complex  (1.6a) 

and 

h(x,t)  =  (hj+hjX+hjt )exp [ia(x-HBt)l  (1.6b) 

The  latter  "galloping"  wall  shape  is  a  surface  wave  which  induces  disturbances 
that  propagate  at  the  freestream  speed.  0=11^  ^  ginn80ldal  „an  sJlape 

(1.6a)  with  a  real  cannot  induce  disturbances  propagating  at  the  freestream 
speed . 


Flow  in  a  wind  tunnel  with  large-amplitude .  non-sinusoidal  wavy  walls  was 
analyzed  in  Ref.  21  also.  For  the  cases  analyzed,  velocity  disturbances  of 
the  form 


u.v  =  [f (y) ,0(y)]exp[io(x-ct)] 


(1.7a) 


with  one  x-wavenumber .  a,  were  induced  by  large  wall  displacements  which  were 
expressed  as  the  Fourier  series 

oo 

h(x,t)  =  £2  Ikn^sinOaCx-ct)  +  h^^cosa^fx-ct)]  (1.7b) 

n=0 


with  the  series  of  x-wavenumbers .  an=na,  n=0,l,2,***.  The  velocities  were 
expressed  analytically,  and  the  wall  shapes  and  Fourier  coefficients  were 
obtained  numerically. 

In  a  companion  study  to  the  present  report,  Rogler  (Ref.  22)  analyzed 
flow  over  a  semi- inf inite  plate  with  traveling  surface  waves,  including  both 
sinusoidal  and  the  more  general  "galloping"  waves  of  eqn.(1.6b).  The  combined 
case  of  freestream  disturbances  and  traveling  wall  waviness  was  also  analyzed. 
Nonperiodic  flow  in  a  corner  induced  by  a  sinusoidal  wall  was  analyzed  by 
Rogler  (Ref.  23)  as  shown  in  Figures  2a, b.  That  solntion  also  illustrates  the 
superposition  of  the  Kelvin-Helmholtz  solution  and  a  standing  wave  solntion. 


1 .4  Solutions  of  the  boundary  layer  equations  for  flows  past  wavy  walls 

As  described  earlier,  Dryden  (Refs.  4-6)  obtained  solutions  of  the  steady 
boundary  layer  equations  for  an  sinusoidal  variation  of  the  freestream 
velocity.  Calculations  were  made  by  a  modification  to  Pohlhausen’s  method. 
Solutions  of  the  boundary  layer  equations  with  sinusoidal  edge  conditions  and 
periodic  pressure  gradient  have  been  obtained  by  Gortler  (Ref.  24) ,  Quick  and 
Schroder  (Ref.  25;  see  also  Fig, A, 10  of  Ref.  6),  Geropp  (Ref.  26),  and  Walz 
(Ref.  27).  Valz  compared  three  solutions  using  approximate  methods  with  the 
finite-difference  solutions  of  Quick  and  Schroder,  where  the  velocity  (^(x) 
was  designated  as  (1)  a  uniform  flow,  (2)  the  sum  of  a  uniform  flow  and  a 
sinusoidal  perturbation,  and  (3)  the  sum  of  a  uniform  flow  plus  a 
nonsinusoidal  perturbation  in  an  intermediate  sector  to  represent  the  effects 
of  a  flaring  section.  Results  without  a  starting  length  were  also  compared. 
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Fannelop  and  Flugge-Lotz  (Ref.  28)  solved  the  compressible  boundary  layer 
equations  for  subsonic  and  supersonic  flow  over  a  wavy  wall.  They  used  an 
implicit  finite-difference  solution  and  addressed  the  problem  of  a  flat  plate 
joined  to  a  wavy  wall.  The  effects  of  boundary  layer  displacement  were 
incorporated  by  assuming  the  sinusoidal  edge  conditions  to  be  those 
corresponding  to  the  combined  effects  of  surface  waviness  and  displacement 
thickness.  The  boundary  layer  flow  was  found,  the  local  displacement 
thickness  was  calculated,  and  the  wall  shape  was  found  Indirectly.  Walz 
(Ref.  27)  also  compares  two  approximate  methods  at  supersonic  Mach  numbers 
with  and  without  heat  transfer.  Interestingly,  Walz  (page  242)  plots  the 
amplitude  e  of  the  sinusoidal  velocity  which  leads  to  laminar  separation 
within  the  first  wavelength.  This  amplitude  varies  from  2.13%  at  ^=0  to 
about  0.5%  at  Mra=4,0,  but  also  depends  on  heat  transfer.  Presumably,  based  on 
Dryden's  calculations,  this  amplitude  would  also  depend  on  the  phase  6  of  the 
wavlness.  The  smallness  of  these  amplitudes  emphasizes  the  sensitivity  of  the 
boundary  layer  to  waviness,  at  least  when  the  mean  pressure  gradient  vanishes, 
**«/ dx=0 , 

Rogers  (Ref.  29)  developed  an  approximate  method  of  accounting  for 
boundary  layer  effects  on  pressure  drag  produced  by  a  compressible  flow  over  a 
wavy  wall. 

Numerical  solutions  for  supersonic  flows  over  single  and  multiple 
sinusoidal  protuberances  were  obtained  by  Polak,  Werle,  Vatsa,  and  Bertke 
(Ref.  30).  Local  separations  and  r eattachments  were  calculated  using  boundary 
layer  equations  where  the  interaction  of  the  boundary  layer  with  the 

isentropic  supersonic  inviscid  flow  was  modeled  by  linking  the  pressure 
gradient  parameter  to  the  local  inclination  of  the  total  displacement  body. 
This  displacement  body  included  the  effects  of  the  wavy  surface  and  the 
boundary  layer  displacement  thickness. 

In  these  references,  the  pressure  at  the  boundary  layer  edge  was 

impressed  on  the  layer,  in  contrast  to  the  studies  of  Section  1.5  where  the 
pressure  is  not  required  to  he  constant  across  the  boundary  layer. 


1,5  Incompressible  and  compressible  parallel-flow  boundary  layers  with 
smal 1-amp litnde  wavv  walls 

Effects  of  the  mean  boundary  layer  along  the  wavy  wall  have  been  studied 
analytically  and  numerically  by  assuming  a  parallel-flow  boundary  layer  and 
seeking  solutions  of  linear  ordinary-differential  eqnation(s)  subject  to 
linearized  boundary  conditions.  These  studies  assume  that  the  disturbance 
flow  is  of  the  form 

u(x,y)  =  f(y)exp(iax) 

v(x,y)  =  0(y)exp(iax)  (1.8a,b) 

and  similarly  for  the  other  disturbance  variables  in  a  coordinate  system 
moving  with  the  surface  wave.  This  solution  form  has  the  same  x-wavenumber  as 
the  surface  waviness,  and  is  compatible  with  the  linearized  boundary 
conditions  along  y=0 .  However,  it  excludes  steady  waves  of  form  (1.3a,b) 
which  may  be  present. 
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Miles  (Ref,  31)  solved  the  inviscid  Rayleigh  equation  for  flow  over  a 
traveling  wave  in  an  analysis  focused  on  predicting  the  energy  transfer 
between  a  parallel— flow  (air)  boundary  layer  and  gravity  wavea  in  water.  The 
rate  at  which  energy  was  tranferred  was  proportional  to  the  profile  curvature 
at  the  critical  point. 

Benjamin  (Ref.  32)  analytically  studied  flow  over  a  stationary  wavy  wall, 
along  with  other  cases  of  traveling  surface  waves.  With  the  critical  point  at 
the  wall,  viscosity  was  retained  in  a  thin  layer  adjacent  to  the  wall  where 
the  mean  velocity  was  assumed  to  vary  linearly.  Miles  (Ref.  33)  also  Included 
viscosity  in  his  analysis. 

Numerical  solutions  of  the  incompressible,  viscous  Orr-Sommerfeld 
equation  with  a  parallel-flow  Blasius  boundary  layer  and  with  a  wavy  wall  were 
obtained  by  Lessen  and  Qangwani  (Ref.  34)  and  Aldoss  and  Reshotko  (Ref.  35). 
Boundary-value  problems  related  to  Refs.  32,  34,  and  35  are  described  in 
Section  1.6  below. 

Inger  and  Williams  (Ref.  36)  solved  the  compressible,  parallel— flow 
equations  for  flow  over  a  stationary  sinusoidal  wall.  The  Prandtl  number  was 
assumed  to  be  unity.  The  inviscid,  nonconducting  solutions  were  matched  to 
solutions  in  the  sublayer  which  accounted  for  viscosity  and  heat  conduction. 
McClure  (Ref.  37)  also  obtained  solutions  for  compressible  disturbances  in 
pseudo-laminar  boundary  layers  with  turbulent  mean  profiles.  He  developed 
some  justification  for  studying  turbulent  boundary  layers  over  wavy  walls  by 
using  the  turbulent  mean  profiles  in  the  linearized,  parallel-flow  equations, 
and  by  using  the  ordinary  viscosity  coefficient  rather  than  some  combination 
of  the  ordinary  and  eddy  coefficients.  Numerical  solutions  of  the  6th-order, 
viscous,  heat-conducting  equations  for  flow  over  stationary  wavy  walls  were 
obtained  by  Lekoudis,  Nayfeh,  and  Saric  (Ref.  38)  .  Mean  boundary  layers 
representative  of  laminar  layers  and  turbulent  layers  were  used  In  their 
calculations. 

The  reader  is  referred  to  Section  1.10  for  related  parallel-flow  studies 
of  the  interaction  of  steady  freestream  disturbances  with  a  boundary  layer 
along  a  flat  plate. 

1.6  Boundary— value  problems  with  wavv  walls 

The  subject  of  this  section  is  the  spatial  evolution  of  disturbances  in 
parallel-flow  boundary  layers  along  wavy  walls.  The  studies  surveyed  in 
Section  1.5  assumed  solutions  of  the  form  (1.8a,b)  which  have  the  same 
x-wavennmber  as  the  surface  waviness.  As  discussed  in  Section  1.1,  however, 
these  solutions  are  not  mathematically  complete.  Any  other  solution  which 
satisfies  the  same  equation,  except  with  homogeneous  boundary  conditions  along 
y  -  0,  can  be  superposed.  There  are  other  steady  solutions,  as  well  as  other 
unsteady  solutions. 

None  of  the  boundary-value  problems  studied  by  Tsugd  and  Rogler 
(Refs,  39a, 39b, 40)  and  Aldoss  and  Reshotko  (Ref.  41)  include  all  of  the 
flat-plate  solutions  and  the  wavy-wall  solution  as  well.  Based  on  the 
linearity  of  the  equation  and  boundary  conditions,  however,  we  can  state  that 
a  complete  set  includes 
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(1)  a  finite  set  ofeigenmodes,  including  the 

Tollmien-Schlichting  stability  wave 

(2)  decaying  standing  waves 

(3)  growing  standing  waves 

(4)  downstream  propagating  continuous  spectrum 

(5)  upstream  propagating  continuous  spectrum 

(6)  wavy  wall  solution  corresponding  to 

Refs.  32,  34  and  35 

Refs.  39a,  39b,  40,  and  41  treated  the  problem  by  Laplace  transform  methods, 
and  formally  decomposed  the  initial  conditions  along  the  y-axis  into  the 
various  waveforms.  Aldoss  and  Reshotko  (Ref.  41)  demonstrated  that  the  effect 
of  the  small-amplitude  wavy  wall  was  to  introduce  two  additional  mathematical 
poles  in  the  Laplace  plane.  These  poles  represent  the  stationary  wavy-wall 
solutions  of  the  Orr-Sommerf eld  equation  calculated  by  Benjamin  (Ref.  32), 
Lessen  and  Gangwanl  (Ref.  34)  and  Aldoss  and  Reshotko  (Ref.  35). 

A  difficulty  with  these  boundary-value  problems  Is  that  the  wavlness  will 
influence  the  boundary  conditions  on  the  steady  fluctuations  specified  along 
the  y— axis.  Additional  information  is  required  to  specify  those  conditions  in 
an  elliptic  problem. 

The  present  work  addresses  the  Issue  of  proper  boundary  condition  when  a 
leading  edge  is  present.  While  each  of  the  solutions  (1-6)  are  independent 
and  superposable,  we  shall  show  that  there  is  a  linking  together  between  the 
steady  waves  (2  and  6)  if  a  leading  edge  is  present.  Instead  of  the 
quarter-plane  problem  ,  Qf  Refs.  39-41,  we  consider  a  whole— plane 

problem  (-»<x<-K»,  -«<y<+oo) .  It  may  he  said  that  a  semi-infinite  wavv  wall 

"excite sr>  several  families  of  waves  with  interrelated  amplitudes ,  phases ,  and 
wavenumbers.  This  process  is  closely  related  to  the  excitation  of  standing 
waves  by  freestream  disturbances,  as  analyzed  by  Rogler  and  Reshotko  (Ref.  9). 

1.7  Flow  over  js  nons inusoldal  bump  treated  as  a,  superposition  of  linear 
solutions  for  flows  over  sinusoidal  walls 

If  the  equations  and  boundary  conditions  for  flow  past  a  sinusoidal  wall 
are  linear,  then  a  superposition  of  the  solutions  is  possible,  of  course. 
Periodic  wall  shapes  can  be  represented  as  Fourier  series,  and  nonperiodic 
wall  shapes  can  be  represented  as  Fourier  integrals.  For  example,  Benjamin 
(Ref.  32)  treated  the  case  of  flow  over  a  bump  as  a  superposition  of  his 
analytical  solutions  of  the  Orr-Sommerfeld  equation  for  flow  past  a  sinusoidal 
wall . 


Generally,  the  flow  over  finite-length  and  semi-infinite  length  plates 
with  wavy  surfaces  cannot  be  represented  as  merely  the  superposition  of  the 
solutions  obtained  with  a  plate  that  extends  from  -»  to  +“>.  As  noted  earlier, 
the  reason  is  that  these  wavy  solutions  do  not  form  a  complete  mathematical 
set,  and  other  waveforms  arise  which  need  to  be  included.  The  present 
analysis  describes  some  of  these  additional  waves  and  obtains  their  amplitudes 
and  phases.  We  note  that  there  are  special  cases  where  the  solutions  for  a 
doubly- inf  ini te  plate  can  be  used  to  construct  the  solution  for  a 
semi-infinite  or  finite-length  plate.  An  example  is  the  symmetric  case  of 
inviscid,  irrotational  flow  over  a  plate  with  waviness  on  the  top  side  which 
is  180  out  of  phase  with  the  waviness  on  the  bottom  side. 
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1*8  Nome rjLcajl  solutions  of  partial  differential  sanations  for  flows  past  wavy 
walls 

Markatos  {Ref.  42)  obtained  numerical  solutions  of  a  set  of  elliptic 
partial  differential  equations  with  constant  density  and  with  periodic 
boundary  conditions  on  the  upstream  and  downstream  boundaries  of  the 
computational  domain 

vU+p,y)  =  v(x,y) ,  u(x+p)=n(x,y)  (1.9a) 

In  a  curvilinear,  orthogonal  coordinate  system,  solutions  are  obtained  by  a 
finite  difference  method.  Flows  were  calculated  with  a  nous inusoidal  surface 
waviness  of  the  form 

h(x)=  -(27i)-l/lexp(-xJ/2)  [C0-C1r+Ca(xl-1)] 

This  surface  represents  water  waves  in  a  coordinate  system  moving  at  the  wave 
speed.  Markatos  considered  both  laminar  and  turbulent  cases,  and  examined  the 
temperature  and  concentration  fields  as  well  as  the  flowfields  in  an 
exceptionally  thorough  study.  The  mean  flows  were  planar,  and  a  two— equation 
turbulence  model  was  used,  where  two  additional  differential  equations  were 
solved  for  the  kinetic  energy  of  turbulence  and  its  rate  of  dissipation. 

Using  Fourier  and  Chebyshev  series,  Balasubramian  and  Orszag  (Ref.  43) 
calculated  steady,  planar  flows  over  sinusoidal  and  non-sinusoidal  periodic 
walls  with  periodic  inflow  and  outflow  boundary  conditions.  The  walls  were 
stationary,  and  both  laminar  and  turbulent  flows  were  studied.  A  one-equation 
eddy  viscosity  model  was  used  in  the  simulations  of  turbulent  flow.  Conformal 
mapping  was  used  to  transform  the  region  above  the  plate  into  a  rectangular 
strip,  then  this  strip  was  mapped  into  a  finite  rectangnlar  domain. 

For  low  Reynolds  number  flows.  Chin  (Refs.  44,45)  developed  and  applied 
conformal  mapping  techniques  to  flows  over  wavy  walls.  Although  we  are  not 
directly  concerned  with  lubrication  problems,  we  cite  these  references  because 
the  mapping  techniques  are  useful. 


1-9  Experiments  with  flows  past  wavy  walls 

Stanton,  Marshall  and  Houghton  (Ref.  46)  experimentally  measured  flow 
oyer  a  wave  train  where  the  amplitude  and  wavelength  grew  in  the  downstream 
direction.  Motzfeld  (Ref.  47)  investigated  flow  over  b  rigid  wavy  surface 
composed  of  three  wavelengths. 

A  wavy  wall  with  waves  propagating  either  downstream  or  upstream  was  used 
by  Kendall  (Ref.  48)  to  study  their  influence  on  the  structure  of  a  turbulent 
boundary  layer  developing  along  that  wall.  Kachanov,  Kozlov,  Kotjolkin, 
Levchenko,  and  Rndnitsky  (Ref.  49)  experimentally  measured  the  laminar 
velocity  profiles  over  stationary  waviness  and  compared  the  results  with 
Gortler's  solutions  of  the  boundary  layer  equations  (Ref.  24).  As  noted 
McClnre  (Ref,  37)  and  Inger  and  Williams  (Ref.  36)  experimentally  and 
analytically  studied  the  compressible  turbulent  flows  above  smal 1-ampl itnde 
stationary  wavy  walls.  Bertram,  et  al  (Ref.  50)  investigated  hypersonic  flows 
past  wavy  surfaces. 
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A  wavy-wall  wind  tunnel  was  built  at  Cambridge  University  and  nsed  for 
studies  with  unsteady  airfoils  and  cascades  (Refs.  51-54).  In  that  facility 
over  a  length  of  2.7m,  a  system  of  cams  and  springs  deformed  flexible  metal 
sheets  to  produce  a  sinusoidal  wall  shape  of  wavelength  1.8o  over  a  central 
region  of  1.2m.  The  remaining  1.5m  were  used  for  the  two  end  regions  where 
the  amplitudes  diminished  to  zero  so  that  the  wavy  wall  would  flare  smoothly 
into  the  stationary  walls  of  the  tunnel. 

Holmes  (Ref.  51)  found  that  significant  discrepancies  existed  between  the 
theory  and  the  experimental  data  when  the  tunnel  operated  in  the  out-of-phase 
mode  where  purely  longitudinal  disturbances  were  introduced  along  the 
centerl ine . 

We  believe  that  these  discrepancies  between  theory  and  experiment  are 
related  to  the  presence  of  standing  waves  we  find  in  the  present  analysis  and 
arise  because  the  sinusoidal  wavy  walls  do  not  extend  forever  upstream  and 
downstream.  In  addition  to  the  traveling  waves  analyzed  by  Horlock  (Ref.  20) 
and  Rogler  (Ref.  21),  one  would  expect  to  find  a  superposition  of  unsteady 
exponentially-decaying  and  growing  standing  waves  of  the  form  (1.3a,b).  A 
neutral  standing  wave  of  form  u(t)=u„exp (-iut ) ,  v=0  could  also  arise  for 
certain  wall  geometries  and  phase  relations  between  the  walls  in  some  sections 
of  the  tunnel. 


1.10  Boundary  layer  over  a  flat  plate  with  a.  steady  sinusoidal  exterior 
boundary  condition 

While  wavy  walls  introduce  disturbances  at  the  wall  and  the  freestream  is 
ordinarily  assumed  to  he  quiescent  in  analytical  studies,  a  related  problem  is 
the  introduction  of  disturbances  from  the  freestream  and  the  plate  is  assumed 
to  be  flat.  If  the  x-wavenumbers  and  phase  speeds  of  disturbances  in  the  two 
cases  are  the  same,  then  many  features  are  quite  similar.  They  have  the  same 
critical  height  y  where  U(yc)-c,  and  they  can  have  similar  viscons  sublayers 
at  the  wall.  Based  on  first-order  boundary  layer  theory,  the  edge  conditions 
and  the  pressure  gradients  may  be  the  same  or  similar.  If  the  disturbances 
are  weak  and  linearized  equations  and  boundary  conditions  are  adequate,  then 
the  two  problems  may  share  certain  superposable  parts  of  their  solutions.  If 
the  disturbances  are  stationary  (c=0) ,  then  the  critical  layer  and  viscous 
sublayer  merge  into  a  single  layer  at  the  wall. 

The  problem  of  a  parallel-flow  boundary  layer  along  a  pl&te.  with 
periodic  fluctuations  at  the  boundary  layer  edge  is  the  steady  limit  (u>  =  0) 
of  the  problem  analyzed  by  Rogler  (Ref.  55) ■  The  boundary  layer  was 
represented  as  a  parallel-flow  Blasius  layer,  and  the  Orr-Sommerf eld  equation 
was  solved  subject  to  the  exterior  behavior 

0(y)  =  Ae+®y  +  Be+yy  +  Ce-ay  +  De_yy 

and  subject  to  the  usual  (flat)  wall  boundary  conditions  of  impermeability  and 
no-slip  at  y-0 .  The  constants  A  and  B  were  specified,  and  the  constants  C  and 
D  were  found  numerically.  This  problem  is  related  to  the  parallel— flow 
treatment  of  incompressible  flow  over  a  wavy  wall  (Refs.  32,34,35). 
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Lighthill  (Ref ■  56)  analyzed  the  inviscid,  nonconducting  case  of  a 
stationary,  weak  distorbanco  impinging  on  a  parallel-flow  boundary  layer  along 
a  flat  plate  with  a  supersonic  freestream.  This  is  related  to  the  cases  of 
compressible  flow  over  a  stationary  wavy  wall  (Refs.  36-38). 

Lighthill  (Ref.  57)  also  considered  the  above  case  except  he  included  the 
effects  of  viscosity  in  a  sublayer  near  the  wall.  The  low  Mach  number  mean 
flow  near  the  wall  was  represented  as  an  incompressible  constant  property  flow 
with  disturbances  governed  by  the  Orr-Sommerf eld  equation  with  a  linear  mean 
profile  and  phase  speed  zero.  Neglecting  terras  of  the  order  of  a* ,  the 
solution  indicated  that  "the  effect  of  the  sublayer  on  the  flow  outside  it  is 
exactly  as  if  the  inner  part  0<y<0.78L  of  it  were  replaced  by  a  solid  wall  and 
the  remainder  were  not  subject  to  viscosity",  where  the  viscous  length  is 

1  and  is  the  wall  value  of  the  kinematic  viscosity.  The 
case  where  the  a  terms  were  included  in  the  analysis  was  also  considered. 

Other  analyses  with  freestream  disturbances  interacting  with 
semi-infinite  plates  which  are  closely  related  mathematically  to  the  present 
work  are  the  interaction  of  an  array  of  square  vortices  with  a  semi-infinite 
plate  (Ref.  7),  seven  other  disturbance  forms  interacting  with  plates 
(Refs.  8),  and  the  representation  of  those  interactions  as  the  superposition 
of  traveling  and  standing  waves  (Ref.  9). 


1 .11  Plan  of  this  stndv 

The  solution  for  incompressible,  irrotational  flow  over  a  semi-infinite 
sinnsoidal  wall  is  obtained  by  conformal  mapping  and  integral  transforms  in 
Section  2.0,  In  Section  3.0,  this  solution  downstream  of  the  leading  edge  is 
represented  as  the  sum  of  the  periodic  Kelv in-Helmholtz  solution  and  a 
secondary  flow  which  arises  because  of  the  leading  edge. 

The  normal  velocity  along  a  line  x<=constant  is  fonnd  above  the  plate,  and  the 
velocity  associated  with  the  Kelvin-Helmholtz  solution  is  subtracted  out.  The 
remaining  velocity  is  Fourier  analyzed  in  Section  4.0  and  represented  as  a 
Fourier  integral  of  waves  of  form  exp(-0x)alnpy.  Upstream  of  the  leading 
edge,  the  flow  is  represented  as  a  superposition  of  growing  waves  of  form 
exp(+0x)cospy  in  Section  5.0  The  study  is  summarized  in  Section  6.0. 
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2  0  FORMULATION  AND  SOLUTION 

By  conformal  mapping,  we  will  obtain  the  flowfield  with  stationary 
surface  waviness  on  a  semi-infinite  plate.  The  method  of  solntion  is  similar 
to  the  procedure  used  to  analyze  the  interaction  of  fieestream  disturbances 
with  a  semi- inf ini te  plate  (Refs.  7-9).  We  are  seeking  a  solution  which 
satisfies  Laplace's  equation,  =  0,  where  J*'*'  is  the  streamf unction 

associated  with  the  flow  induced  by  the  surface  waviness.  For  sign 
convention,  the  streamfunction  is  related  to  the  longitudinal  and  normal 
velocities  by  -5*"*'=  u^  and  respectively. 

The  geometry  and  notation  for  the  semi-infinite  plate  with  a  wavy  surface 
is  shown  in  Figure  3.  The  displacement  of  the  wall  is  the  shifted  sine  wave 

h(x,0)  =  h#sin(ax-8)  (2,1) 

where  h0  is  the  amplitude,  0  is  the  phase  angle  of  the  waviness,  and  a  is  the 
x-wavennmber . 

The  linearized  relation  between  the  surface  displacement  and  the  normal 
velocity  along  the  x-axis  is 

vU)  <x20,0)  =  U^dh/ax  (2.2) 

-  U^aiijjcos  (ax-0)  =  (for  y=0,x20)  (2.3a) 


The  streamfunction  is 

^  (x20,0)  =  Ugghg sin (ax-0)  +  (2.3b) 

where  the  streamfunction  which  is  a  constant  or  function  of  time,  can  be 
neglected  with  this  semi-infinite  plate.  These  boundary  conditions  are  valid 
only  downstream  of  the  leading  edge. 

Hereafter,  the  velocities  are  nondimensional ized  against  o^h,, ,  the 
streamfunction  is  nondimensional  ized  against  JiDw>h0 ,  and  the  x  and  y 
coordinates  have  been  nondimensional ized  against  the  half-wavelength  k/2  of 
the  sinusoidal  surface  wavlness. 

For  the  case  of  the  waviness  being  in-phase  on  opposite  sides  of  the 
plate,  the  streamfunction  will  be  symmetric  about  the  x-axis.  Hence,  the 
homogeneous  Neumann  boundary  condition  is  applicable  npstreao  of  the  plate 

=  0  (for  y=0,  x<0)  (2.4) 


The  solution  is  assumed  to  vanish  far-upstream  from  the  plate  and  at 
great  distances  laterally  from  the  plate.  Far-downstream  of  the  leading  edge, 
one  would  expect  that  the  Kelvin-Helmholtz  solution  would  be  recovered. 
Hence,  we  require  tbat  the  solution  be  bounded  as  jr9+®>.  These  exterior 
boundary  conditions,  along  with  the  Dirichlet  and  Neumann  conditions  along  the 
x— axis  and  Laplace's  equation  yield  a  well— posed  elliptic  mathematical  system. 
The  streamfunction  will  be  symmetric  about  the  x-axis,  so  we  can  confine 
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attention  to  either  half-plane  which  is  above  or  below  the  plate. 

We  conformally  map  the  half-plane  above  the  x-axis  onto  the  quarter-plane 
as  shown  in  Figure  4  by  the  mapping,  p  =  zl/j,  where  z  =  x+iy  and  p  =  5+i,,. 
The  conformal  mapping  alleviates  the  problem  of  split  boundary  conditions 
along  the  x-axis,  with  a  Neumann  boundary  condition  for  x<0  and  a  Dirichlet 
boundary  condition  for  x>0  (Eqns.  2. 3b, 2. 4).  Under  this  mapping,  the 
resultant  mathematical  system  is  defined  in  Figure  4. 

The  solution  is  obtained  by  a  cosine  integral  transform  in  the  £ 
direction.  The  details  are  presented  in  Ref.  7.  The  solution  for  the 
streamf unction  of  the  flaw  induced  by  the  semi-infinite  plate  with  a 
sinusoidal  wall  is 

=  tcos0  Real[-oosjizS4(-nz)-sinitzC2(-7Tz) 

+cosrrzC2  (-nz)  -sinnzSa  (— jtz)  +sinnz] 

-sin©  Real [-cosnzS, (-nz)-sinnzC, (-nz) 
~C0snzCj{-CTz)+sxnnzS2(-Tiz)+C0  9JtzJ}/i!  (2.5) 

where  S2  and  C2  are  the  Fresnel  integrals 


Sj(z)  =  (2n)  1^*Jt  sin ( t )  dt  (2.6a) 


c2(z)  =  (2n)  1^*Jt_:lJ,1cos(t)  dt  (2.6b) 


The  Fresnel  integrals  are  evaluated  numerically,  either  as  asymptotic  series 
valid  for  large  argument  or  series  expansions  valid  for  small  argument 
(Re f  *  5 8)  , 

The  velocities  corresponding  to  solution  (2.5)  are 

u<i)  =  -^U)  =  _iF(2jQ) 

v(i)  =  p(i>  =  F(Z,Q) 

where  F  is  the  complex  function 

F  =  sin<Ttz-e)[S2(-jiz)-C2(-jiz)] 

— cos (nz-0) [SJ(-nz)+C2(-nz)-l] 

-(sin6+cos0)/[n(-2z) 1/1] )  (2,7c) 

For  large  |z|,  computational  experience  has  shown  that  it  is  better  to  recast 
the  solution  as 


(2.7a) 

(2.7b) 
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Edge 


Figure  3.  Geometry  and  notation  for  flow  past  a  wavy,  semi-infinite  plate. 


=— -sin(jr£  -  ©)  + 


Figure  4.  Conformal  mapping  of  the  upper  half-plane  onto  the  quarter-plane. 
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F  =  sin©[f  -  g  -  1/ [ir<-2z)  ]  J 

+  cos6[f  +  g  -  1/[k(-2z)i/j]]  (2 ,7d) 

where  f(z)  and  g(z)  are  related  to  the  Fresnel  integrals  by 

f(z)  =  -  S( z) ] cos (nz* /2)  -  —  C(  z)  ]  sin  (rrz  4/2) 

g(z)  =  C(z)]cos<7tz1/2)  +  SCz>]sin(iizV2) 

S(z)  =  S1{„z1/2)i  C(z)  =  Cj (nz*/2) 

The  asymptotic  series  for  f  and  g  (Ref*  58,  page  302) 

00 

?Tzf(z)~  1  +]H(-l>m  [l-3---{4nr-l)/(nzJ)anl 
o=i 

00 

nzg( z)  ^(-1)“  [1*3*  •  •  ( 4m+l) / (jizs ) 

iifo 


converge  rapidly.  The  solntion  form  (2.7d)  avoids  the  problems  of  calculating 
the  small  differences  between  large  numbers.  While  the  direct  evaluation  of 
solution  (2.7b)  is  limited  to  abont|z|=1.8  when  32  bit  arithmetic  is  used,  and 
is  limited  to  abont|z|=5.2  when  64  bit  arithmetic  is  used,  solution  (2.7d)  has 
been  evaluated  at|zt=400,  which  is  200  wavelengths  from  the  leading  edge. 

Figures  5a, b  show  the  disturbance  streamline  pattern  for  the  flow  induced 
by  the  semi- inf inite  wavy  wall  with  two  phase  angles.  These  figures,  and  the 
others,  do  not  include  the  uniform  mean  flow.  The  induced  flow  depends  on  the 
phase,  O,  of  the  surface  wave.  The  wavy  surfaces  on  opposite  sides  of  the 
plate  are  in  phase.  One  can  clearly  see  the  upstream  influence  of  the 
waviness,  as  expected  since  Laplace's  equation  is  elliptic. 

A  singularity  occurs  at  the  leading  edge  for  the  velocities,  with  the 
validity  of  the  theory  consistent  with  linearized  airfoil  theory.  The  region 
where  the  velocities  are  large  can  be  made  as  small  as  desired  by  making  the 
wall  waviness  smaller.  Hence,  this  is  the  proper  first— ordex  solution  as  long 
as  the  flow  is  onseparated  at  the  leading  edge.  Since  we  are  interested  in 
the  practical  case  of  flows  past  bodies  with  streamlined  leading  edges,  the 
assumption  of  onseparated  flow  implicitly  contained  in  the  Laplace’s  equation 
appropriately  models  that  feature  without  unduly  complicating  the  analysis 
with  details  of  the  nose  geometry  and  plate  thickness. 

Far-downstream  of  the  leading  edge,  it  appears  in  Figures  5a, b  that  the 
tel vin-Helmhol tz  solution  is  recovered.  However,  the  Kel vin-Helmholtz 
solntion  is  one  part  of  the  solution  for  all  xiO+.  We  shall  demonstrate  this 
feature  in  the  next  section  by  rearranging  the  solntion. 
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(5a)  6  =  0 


(5b)  e  =  n/2 


Figure  5.  Streamlines  for  the  flow  induced  by  surface  waviness  on  a 
semi-infinite  plate  for  two  phase  angles  of  the  surface  waviness. 
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3,0  A  SEPARATION'  OF  THE  SOLUTION  DOWNSTREAM  OF  THE  LEADING  EDGE 

INTO  TWO  CLASSES  OF  WAVES 

Recent  analyses  (Refs.  39a, 39b, 40, 41)  which  have  identified  the  basic 
solutions  of  the  Orr— Sommerf eld  equation  offer  guidance  on  how  to  decompose 
solution  (2. 5, 2. 7a, b)  into  its  basic  Fourier-Laplace  constituents.  Rather 
than  resort  to  a  formal  boundary-value  problem  in  the  quarter-plane  x20,y20 
for  this  inviscid  problem  with  a  uniform  mean  flow,  we  will  list  the 
simplified  versions  of  the  waves  discussed  in  Section  1.6.  The  eigenmodes  and 
the  upstream-propagating  continuous  spectrum  do  not  appear  if  the  mean  flow  is 
uniform  and  viscous  effects  are  neglected.  For  the  region  downstream  of  the 
leading  edge,  the  Rayleigh  equation  with  a  uniform  mean  flow  is 

(1  --£)(D*-fca)0(y)  =  0  (3.1) 

where  the  normal  velocity  and  its  amplitude  are  related  by 
v(x,y, t)=0(y)exp (ikx-iwt ) .  The  x-wavenumber  can  be  complex.  Solutions  are 
assumed  to  be  bounded  as  y  ■>  ».  The  boundary  condition  along  y=G  depends  on 
the  wavenumber  and  frequency,  and  we  shall  note  the  appropriate  condition  for 
each  wave.  The  formal  boundary-value  problems  properly  sort  out  these 
solutions  and  conditions  through  the  use  of  Fourier  and  Laplace  transforms, 
but  for  brevity,  a  more  direct  statement  will  be  made  here  for  the  subject 
problem. 

The  Rayleigh  equation  (3.1)  only  admits  solutions  of  four  forms  which  are 
bounded  as  y9at>: 

(1)  The  (steady)  classical  Relvin-Helmholtz  fluctuating  flow  with  phase 
speed  zero.  Only  this  solution  has  a  nonhomoge neons  boundary 
condition  along  the  x-axis,  and  the  wavenumber  k  Is  the  same 
wavenumber  as  the  sinusoidal  wall. 

(2)  Irrotational ,  exponentially— decaying  standing  waves,  of  the  form 

v  =  sin(Py)exp(-px) 

where  p  is  real  and  the  wavenumber  is  pure  imaginary,  k=ip.  This 
solution  satisfies  the  usual  impermeability  condition  along  a  flat  plate. 
At  this  point,  we  do  not  know  the  value  of  p.  Later,  we  shall  show  that 
there  is  a  spectrum  of  these  waves. 

(3)  Irrotational,  exponentially-growing  standing  waves  of  the  form 

v  =  sin (py) exp (+px) 

These  growing  waves  must  be  excluded  downstream  of  the  leading  edge  of 
this  semi-infinite  plate,  unless  some  influence  originating  from  positive 
infinity  is  to  be  included.  Growing  waves  of  form  v(x,y)=cos(Py)erp(+px) 
do  exist  upstream  of  the  plate,  however,  as  discussed  in  Section  50. 

(4)  Vortical  fluctuations  convecting  downstream,  which  can  be  ruled  out  here 
since  no  vortical  freestream  disturbances  have  been  introduced  in  this 
problem,  and  none  can  diffuse  from  the  wall  in  this  inviscid  analysis. 
These  fluctuations  would  satisfy  impermeability  on  a  flat  plate. 


22 


AEDC-TR-33- 10 


Hence,  for  *>0,  we  write  solution  (2.5)  as  the  sum  of  two  parts 
i)  =  ^(kh)  +  ^;(s  ) 

where  j^(kh)  -  ^sin(1IX-0)e-,ty  (Kelvin- Helmholtz  solution)  (3.2a) 

*(•)  =  -  ^(kh)  (Secondary  flow)  (3.2b) 

The  secondary  flow  represents  the  influence  of  the  leading  edge  of  the  plate. 
Hiis  secondary  flowfield  satisfies  the  impermeability  condition  for  a  flat 
Plate.* 


The  streamfunction  has  been  evaluated  downstream  of  the  leading 

edge,  and  the  Kelvin-Helmholtz  solution  has  been  subtracted  out.  The 
streamlines  of  the  secondary  flow  which  remains  are  plotted  in  Figures  6a-h 
for  eight  different  phase  angles.  Note  in  Figures  6d-6g  that  the  flow 
reverses  itself  near  the  wall,  and  that  the  size  of  the  cell  along  the  x-axis 
depends  on  the  phase  angle  of  the  surface  waviness.  The  longitudinal  velocity 
along  the  z-azis  is  the  function  u^s^(x,G)  which  appeared  in  eqn.  1.5. 

The  normal  velocity  associated  with  this  second  flow,  v^s^  =  <^s^along 
the  line,  xx  =  0.5,  is  plotted  in  Figure  7  for  two  phase  angles.  These 

velocity  profiles  were  calculated  along  the  dashed  lines  drawn  in  Figures 
6a, f.  The  maximum  y-value  is  much  larger  than  for  the  other  figures  so  that 
the  decay  can  be  illustrated.  Note  that  this  second  flow  satisfies 
impermeability  (v^s^=0  at  y=0)  in  contrast  to  the  nonzero  value  for  the 
Kelvin-Helmholtz  solution. 


4.0  THE  SPECTRUM  OF  STANDING  WAVES 


It  is  often  useful  to  express  the  solution  in  terms  of  its 
Four ier-Lapl ace  constituents,  each  which  can  be  calculated  and  the  waves  then 
superimposed.  The  Kelvin-Helmholtz  solution  (3.2a)  is  already  in  the  form  of 
a  neutrally  oscillating  wave.  The  normal  velocity  for  the  secondary  flow  will 
be  decomposed  by  taking  the  half-range  sine  transform 


with  inverse 


v(s)(z1,p,«)  =  (2/7t)l/1|v(s)(x1,y,G)sinpy  dy 

00 

^Xj.y,#)  =  (2/7t)l/lJv(s)(x1,pfe)sinPy  dp 


f(s 


<4. la) 


(4.1b) 


Since  each  of  the  waves  decay  exponentially  in  the  x-direction,  the  spectrum 
at  position  x  is  related  to  the  spectrum  at  position  hy  the  relation 

v<s>(p,x,e)  =  v(s) (p,z1,0)exp[-p(z-z1)] 

for  x.ZjX).  Hence,  for  other  x-positions,  the  velocity  can  be  expressed  in 
terms  of  the  spectrum  at  by  the  relation 
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(6c)  0  =  0.30n 


(6d)  0  =  0.35n 


Figure  6.  Streamlines  for  the  standing  wave  pattern,  for  various  phase 

angles  of  the  surface  waviness.  This  flow  represents  the  alteration 
to  the  Kelvin-Delmholtz  solution  produced  by  the  leading  edge.  It 
is  composed  of  a  superposition  of  the  basic  waves  shown  in  Figure  8. 
Other  phase  angles  are  shown  in  Figures  6e-h.  The  incremental  value 
of  the  streamf unction  has  two  values,  with  one  value  ten  times 
smaller  than  the  other  to  better  illustrate  the  flow. 
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( 6 g )  ©  =  0.60n 


Figure  6.  Concluded. 
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Figure  7,  The  normal  velocity  fluctuation  for  the  superposition  of  standing 
waves  for  two  phase  angles.  These  are  the  velocity  fluctuations 
along  the  dashed  lines  drawn  in  Figures  6a  and  6f.  The  maximum 
ordinate  is  much  largor  for  this  figure  than  for  the  other  figures 
to  illustrate  the  slow  decay. 
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v*s^(x,y,0)  = 

Oft 

(2/n)1/*jv^s^(p,i1,€))e~^x-Xl^sin0y  dp  (4.1c) 


The  streamlines  for  exponentially-decaying  standing  waves  are  plotted  in 
Figs.  8a-c  with  three  values  of  the  x-waveumber ,  k  =  ip.  The  three  values  of 
k  are  approximately  the  values  corresponding  to  the  maximum  Fourier  amplitudes 
in  Figure  9, 

The  Fourier  spectrum  of  the  velocity  plotted  in  Figure  7  is  shown  in 
Figure  9.  The  spectrum  illustrates  that  the  velocity  v^s^  is  composed  of  a 
superposition  of  steady  waves,  each  of  form 

s ^ (x.y.p ,0)  =  v^®) (p,©)e“PIsinPy  (4.2) 

This  flow  is  also  an  irrotational  flow.  Unlike  the  classical  Kelvin- Helmholtz 
solution,  however,  this  solution  decays  exponentially  downstream  of  the 
leading  edge,  and  satisfies  the  impermeability  condition  at  the  surface  of  a 
fiat  plate.  Although  this  wave  oscillates  neutrally  as  y^  ®  ,  we  have  shown 
that  such  waves  can  originate  from  the  wall.  The  superposition  of  these 
waves,  however,  yields  a  flow  which  vanishes  far  away  from  the  plate. 


5  0  SOLUTION'  UPSTREAM  OF  THE  LEADING  EDGE 
REPRESENTED  AS  A  SUPERPOSITION  OF  GROWING  STANDING  WAVES 


Solution  (2.5)  for  the  streamfunct ion  j£(i)  solutions  (2.7a,b,c)  for 
the  velocities  are  valid  both  upstream  and  downstream  of  the  leading  edge. 
However,  the  representation  =  v^kJ+yts)^  with  v^s^  expressed  in  equation 
(4.1c)  as  a  superposition  of  decaying  standing  waves,  is  valid  only  downstream 
of  the  leading  edge.  What  is  the  corresponding  representation  upstream  of  the 
leading  edge  in  the  half— plane  x<0,  -oo<y<«>7 


Upstream  of  the  leading  edge,  there  is  no  constraint  that  the  solution  of 
the  Rayleigh  equation  (3.1)  satisfy  the  impermeability  boundary  condition. 
Hence,  the  Kelvin-Helmholtz  solution  does  not  appear  in  the  upstream 
half-plane.  Hence  the  possible  solutions  for  are  growing  and  decaying 
standing  waves  and  vortical  fluctuations  of  the  forms 


Standing  waves 

(a)  exp  (+px)cospy 

(b)  exp(+Px)sinpy 

(c)  exp  (-px) cosjly 

(d)  exp (-px) sinpy 


Rectangular  arrays  of 
vortices  convecting  downstream 

(e)  siny(x-t)sinpy 

(f)  siny(x-t)cospy 

(g)  cosy(x-t)sinpy 

(h)  cosy(x-t)cospy 


(5.1) 


Of  these  possibilities,  the  form  (b)  can  be  eliminated  by  the  homogeneous 
Neumann  condition.  Forms  (c,d)  blow  up  as  x  4  -oo  and  can  be  eliminated.  The 
last  four  solutions  (e-h)  represent  vortical  freestream  disturbances  convected 
downstream,  and  they  do  not  appear  in  the  present  problem.  Only  the  waveform 
(a)  remains.  Hence,  the  normal  velocity  upstream  of  the  leading  edge  can  be 
represented  as  the  cosine  integral  transform 
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«•  Streamlines  for  three  examples  of  the  basic  Fourier-Laplac 
constituents  of  the  flow  illustrated  in  Figures  6a-h.  These  at 
stationary  waves  which  decay  like  sin(py)exp(-0x) .  The  flow  i 
ignres  6a  h  represent  a  superposition  of  these  waves,  wit 
spectrums  plotted  in  Figure  9.  These  three  values  of  th 

aPPro*lni*tely  to  the  local  maximums  of  th 
rourier  amplitudes  in  Figure  9. 
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with  inverse 


oo 

,<s)(*x.P.«)  =  (2/n)l/lJ>v(s)(x1,y,e)cospy  dy 

00 

,(s>(x1.y-0)  =  <2/n)1/l|v(s)(*1,p,0)cos^y  dp 


(5.2a) 


(5.2b) 


Of  coarse  for  an  even  function  in  y,  a  cosine  transform  is  adequate.  The  fact 
that  they  are  growing  standing  waves  appears  in  the  next  step.  These  are  the 
Fourier  amplitudes  and  velocity  along  some  line  x=x1<0.  Since  each  of  the 
waves  grows  exponentially  in  the  streamwise  direction,  the  amplitude  at  x  is 
related  to  the  amplitude  at  x1  by 

v(x,pt0)=v(i1,jj,0)exp  [+p(x-xj)  ] 

Hence,  the  velocity  at  the  position  x,y  is 

v^ s  ^ (x,y,0)  = 
oo 

=  (2/jr)  l/aJv{s)(xl,p,6)e+P(x“xi)cospy  dp  (5.3) 

0 

v 

for  x,  xx<0 .  There  is  no  "wavy— wall"  solution  for  x<0.  Hence,  v^^=v^s^  for 
x<0.  The  cosine  integral  expresses  the  solution  as  a  superposition  of 
exponentially-growing,  symmetric,  standing  waves. 


6.0  SUMMARY,  DISCUSSION  AND  CONCLUSIONS 

Conformal  mapping  of  the  irrotational  flow  Induced  by  surface  waviness  on 
a  semi-infinite  plate  has  yielded  solution  (2.7b)  for  the  normal  velocity 
which  is  valid  both  upstream  and  downstream  of  the  leading  edge.  However, 
downstream  of  the  leading  edge,  it  is  instructive  to  express  this  solution  in 
terms  of  two  families  of  waves: 

v^(x,y,6)  =  cos(jtx-4J)exp  (-ny) 

00 

+  (2/n)l/,|v(s)(x;l,p,e)sinpy  exp[-p(x-x1)]  dp  (6.1) 


(1)  The  classical  Kelvin— Helmholtz  flow  induced  by  a  wavy  wall,  with 
streamlines  plotted  in  Figure  1.  This  solution  satisfies  a 
nonhomogeneons  wall  boundary  condition  along  the  x-axis. 

(2)  A  secondary  flow  which  can  be  represented  as  a  superposition  of 
waves  which  decay  exponentially  in  the  streamwise  direction. 
This  flow  accounts  for  the  presence  of  the  leading  edge  of  the 
semi-infinite  plate.  This  snectrnm  has  resulted  from  surface 
waviness  which  has  only  one  x-wavenumber.  These  waves  and  their 
combined  pattern  satisfy  the  wall  boundary  condition  on  a  flat 
plate,  vCs'=0  at  y=0,  although  they  have  been  excited  by  the 
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combined  effects  of  the  wall  waviness  and  the  leading  edge.  The 
Fourier  amplitudes  v's'(x1(y,©)  are  plotted  in  Figure  9. 
Far-downstream  of  the  leading  edge,  the  secondary  flow  vanishes 
and  the  only  solution  which  remains  is  the  classical 
Kelvin-Helmholtz  solution. 

Also  downstream  of  the  leading  edge,  it  is  instructive  to  write  the 
streamwise  velocity  along  the  plate,  in  dimensional  form,  as 

uoo<*)  =  n.  +  u(kh}Cx,o,e)  +  u(s) Cx,o,e)  (6.2) 

This  expression  emphasizes  that  the  longitudinal  velocity  along  the  plate  is 
composed  of  the  value  it  would  have  if  the  plate  were  flat,  a  sinusoidal 
perturbation  associated  with  the  Kelvin-Helmholtz  solution,  and  a  secondary 
flow  consisting  of  a  superposition  of  exponentially-decaying  waves. 

These  simplified  forms  of  the  waveforms  are  based  on  a  uniform  mean  flow 
which  is  inviscid.  Table  1  summarizes  the  Rayleigh  solutions  (with  U=l)  and 
the  Orr-Sommerf eld  solutions  for  the  two  waveforms  described  above. 

The  flowfield  upstream  of  the  leading  edge  can  be  represented  as  a 
superposition  of  symmetric,  growing  standing  waves 

y(i)  =  V(s)(X,y,e)  = 


-(2/n)  1/,4Jv^s^  cosPyexp  E+lHx-Xj)  ] 


dp 


(6.3) 


Exponentially-decaying  waves  of  the  unsteady  form 

v(x,y,t)  =  0(y)exp(-px-iwt)  (6.4) 

with  0(y)=sinpy  also  have  been  found  previously  in  analyses  of  vortices 
colliding  with  a  semi-infinite  plate  (Ref.  9).  Unsteady  waves  of  the  more 
general  form  (6.4)  have  also  been  found  as  solutions  of  the  Orr-Sommerf eld 
equation  (Ref.  2) .  This  analysis  illustrates  that  steady  versions  of  these 
standing  waves  originate  in  the  flow  abont  a  wavy  semi-infinite  plate. 
Further  study  is  required  to  describe  the  effects  of  viscosity  and  the  mean 
boundary  layer.  While  these  waves  oscillate  neutrally  as  y*®  >  this  analysis 

shows  that  they  can  originate  from  the  wall. 


The  viscous  form  of  the  steady  exponentially-decaying  wave, 
v  =  0(y) exp (— Px) ,  has  not  been  calculated.  The  unsteady  solutions  for  both 
decaying  and  growing  standing  waves  in  boundary  layers  as  solutions  of  the 
Orr-Sommerf el d  equation  are  documented  in  Ref.  2.  When  the  frequency 
vanishes,  however,  the  exterior  solution  takes  on  the  form  of  four  neutrally 
oscillating  waves 


0(y)  =  Ae1^  +  Be_ify  +  Ceipy  +  De_iPT 
where  y=R1^2 (p+p2/Rg) ,  in  contrast  to  the  exterior  solution 
0(y)  =  Ae"“*  4  Ce^y  4  De_iPy 


(6.5a) 

(6.5b) 
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where  m  (  0Rg  p  iwRg)  ^  for  the  unsteady  case.  Numerical  solutions  for  very 
small  frequencies  suggest  that  a  very  slowly  decaying,  high  frequency 
oscillation  appears  and  survives  into  the  freestream.  In  the  exact  steady 
limit,  these  oscillations  may  take  the  form  exp(+iyy)  as  noted  above.  These 
additional  steady,  viscous  solutions  warrant  further  study. 

This  study  has  assumed  a  semi-infinite  plate.  We  believe  that  the 
secondary  flow  generally  will  be  stronger  for  finite-length  plates,  than  for 
semi-infinite  plates  with  wavy  walls,  at  least  for  certain  phase  angles  and 
plate  lengths.  This  conjecture  warrants  detailed  study. 

Finally,  we  recognize  that  there  is  a  relationship  between  the  solutions 
for  the  cases  of 

(1)  freestream  vorticity  disturbances  convecting  downstream  and 
encountering  a  plate 

irrotational,  traveling  freestream  disturbances  propagating  at 
speeds  different  from  0^  and  encountering  a  plate 

flow  past  a  semi-infinite  plate  with  stationary  wall  waviness 


Table  1 

COMPARISON  OF  SOLUTIONS  FOR  FLOW  PAST  WAVY  WALLS 


(a)  Solutions  of  Rayleigh  Equation  (b)  Solutions  of  Orr-Sommerf eld  Equation 
with  a  Uniform  Mean  Flow  with  a  Nonuniform  Mean  Flow 


(2) 

(3) 


(la)  Kelvin-Helmholtz  Solution  (lb) 

varying  as 

v(x.y)=v0exp(-py)sin(px) 
which  satisfies  nonhomoge neons 
boundary  condition  at  y=0, 

(See  Figure  1  and  Ref.  1) 

(2a)  Steady,  Decaying  Solution  (2b) 

varying  as 

v(x.y) =v0sin(Py)exp(-px) 
which  satisfies  homogeneous 
boundary  conditions  at  y=0. 

(See  Figures  8a, b;  Figure  6a— h 
is  a  superposition  of  these 
waves) 


Benjamin  (Ref.  32),  Lessen— Gangwani 
(Ref.  34) , end  Aldoss-Reshotko  (Ref.  35) 
solution  of  form  v(x,y)=0(y)exp(iux) 
which  satisfies  nonhomoge neons 
boundary  condition(s)  at  y=0. 

Limiting  Steady  Case  (w=0)  of  the 
Exponentially— Varying  Standing  Wave 
varying  as 

v(x,y,  t)=0(y)exp  (— px-iint ) 
which  satisfies  homogeneous 
boundary  conditions  at  y=0 . 

(The  unsteady  case  has  been 
calculated  (Ref.  32);  the  steady 
case  has  not  been  calculated) 


Stokes  solution  for  an  unsteady 
surging  of  the  freestream  is  the 
p=0  limit  of  the  unsteady,  standing 
waves 
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(the  case  considered  in  this  study) 

(4)  flow  past  a  semi-inf inite  wavy  wall  with  traveling  surface  waves 

The  relationships  between  these  flows  for  seemingly  rather  different 
cases  of  freestream  disturbances  and  wall  waviness  warrant  further  study  and 
development. 
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English 

c-u>/a 

C(z> 

C*(z> 
f  (z) 

F(z,e) 
g  (z) 

h(x) 

K 

h±  ,h2 

i 

k 

S  ( z ) 

S2(z) 

t 

u,v 


y 

z=x+iy 


NOMENCLATURE 

phase  speed 

z 

Cosine  Fresnel  integral,  C(z)  =  |cos(jtt J/2)dt 

C(z)=C1(nz2/2) 

z 

Cosine  Fresnel  integral(2),  C2  (z)={27t)-:l/I  jt_1^*costdt 
function  related  to  the  Fresnel  integrals  by 

f(z)  =  S(  z)  ]  cos  (itz2/2)  -  [-J--  C(z)lsin(nz*/2) 

function  defined  in  eqn.(2.7c) 
function  related  to  the  Fresnel  integrals  by 

g(z)  -  [-£  C(z)]cos(nz  /2)  +  [-£ —  S(  z)  ]  sin(xtz1/2) 
displacement  of  the  wavy  wall 

amplitude  of  the  sinusoidal  wall  waviness 

amplitudes  of  the  galloping  surface  as  defined  in  eqn.fl.tfb) 
(-1>1/4 

x-wavenumber  in  the  Rayleigh  eqn.O.l) 

Sine  Fresnel  integral,  S( z)  =  |sin(nt*/2)dt 
S( z) “S^ (nz* /2 ) 

Sine  Fresnel  integral{2),  S2 (z)=(2n)-l/aJt-l/1sintdt 
time 

disturbance  velocities  in  the  x  and  y  directions 

mean  x-velocity  in  the  freestresm  that  would  exist 
if  the  plate  wore  flat 

coordinate  parallel  to  plate  and  in  streamwise  direction 
coordinate  normal  to  the  plate 
complex  coordinate 
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Greek  and  Script 
a  x-wavenumber 

P  y-wavenumber 

V2  Laplacian  operator 

e  small  quantity  proportional  to  the  amplitude  of  the  surface 

waviness 

X  wavelength  of  the  surface  waviness 

^  *■)  disturbance  streamf unct  ion  induced  by  the  surface  waviness 

on  a  semi-infinite  plate 

Kelvin— Ilelmholtz  streamfunction  for  flow  past  a  wavy  wall 

streamfunction  composed  of  a  pattern  of  standing  waves 
and  representing  the  alteration  to  caused  by  the 

leading  edge 

p  q+in»  complex  coordinate  in  the  conformally  mapped  plane 

0  phase  angle  of  the  surface  waviness 

0(y)  complex  amplitude  of  the  normal  velocity  fluctuation 

frequency 

4,ti  coordinates  in  the  conformally  mapped  plane 

Superscripts  ,  Subscripts .  and  Miscellaneous  Notation 

i  associated  with  the  fluctuating  flow  induced  by  surface 

waviness  on  a  semi-infinite  plate 

kh  associated  with  the  Kelvin-Helmholtz  solution 

Real[  ]  real  part  of  [  ] 

s  associated  with  the  secondary  flow  which  is  generated  by  the 

combined  effects  of  the  surface  waviness  and  the  leading  edge 

Characteristic  quantities  used  in  the  nondimens ional izat ions 

X/2  characteristic  length 

cth,!^  characteristic  disturbance  velocity 

bgp  characteristic  mean  velocity 

ah0Ueo/n  characteristic  disturbance  streamfunction 
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